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Transport of Brownian particles on a simple sawtooth potential subjected to both unbiased thermal 
and nonequilibrium symmetric three- level Markovian noise is considered. The new effects of three 
and four current reversals as a function of temperature are established in such correlation ratchets. 
The parameter space coordinates of the fixed points associated with these current reversals and the 
necessary and sufficient conditions for the existence of the novel current reversals are found. 
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The aim of this paper is to study current reversals 
in controlled transport of Brownian particles (!]] in- 
duced by symmetric nonequilibrium noise in ratchets 
with a simple sawtooth potential, bearing in mind po- 
tential applications for separation of nanoobjects [||. 
It is known that current reversals in ratchet systems 
can be engendered by changing various system param- 




including the flat- 
the correlation 
, the temperature 



eters UlifflJIBgg 

ness parameter of the noise 

time of nonequilibrium fluctuations 

in multinoise cases [|ll[ , the power spectrum of the noise 
source j jl2| |, the shape of the potential [0|, the number 
of interacting particles per unit cell [|l4f, and the mass 
of the particles fill . As a rule, these results have been 
obtained either at the limits of slow and fast noises or by 
numerical methods. At the same time, analytic results 
would greatly facilitate the study, especially in the inter- 
mediate regimes of the system parameters which is the 
realm of biology. It is especially difficult to obtain ana- 
lytic expressions for the multinoise ratchets. A fortunate 
exception here is the symmetric three-level telegraph pro- 
cess (trichotomous noise) which is rich enough physically 
and can at the same time be treated analytically P,p|,[L6| . 
In this paper we will prove, on the basis of the leading 
order term of a series expansion for the probability cur- 
rent in terms of inverse flatness, that there exist three 
and four current reversals of the probability current as a 
function of temperature D. Never before have more than 
two current reversals with D been reported for correla- 
tion ratchets with a simple sawtooth potential. (At the 
same time, in the case of rocking ratchets, infinitely many 
current reversals may occur |l7|].) We will also derive the 
necessary, and the necessary and sufficient conditions for 
the existence of these novel current reversals. 

A zero-mean trichotomous Markovian stochastic pro- 
cess Z(t) consists of jumps between three values z = 
{a, 0, —a}, a > 0. The jumps follow in time the pattern of 
a Poisson process, the values occuring with the stationary 
probabilities P s (a) = P s (—a) = q and P s (0) = (1 — 2g), 
where < q < 1/2. Denoting the state space and the 
transition matrix of our trichotomous process, respec- 
tively, by {ai}:={a, 0, — a} and Tij := P{a^, t +r |aj, t}, 
z, j = 1, 2, 3, r > 0, we have 
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FIG. 1. Four current reversals vs temperature D for fixed 
d = 0.0045 and a = 27.75. The switching rates are respec- 
tively (1) v x = 649.651630, (2) v 2 = 646.770000, and (3) 
V3 = 643.816291. Curve (2) has four single zeros. Curves (1) 
and (3) have two single zeros and one two- fold zero. (At the 
fixed values of d and a, for curves with v values smaller than 
vz and greater than v\ there can occur at most two current 
reversals.) 
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where v >0. In a stationary state, the fluctuation process 
satisfies (Z(t)) = and (Z(t + r)Z(t)) = 2qa 2 exp(-z/r), 
where the switching rate v is the reciprocal of the noise 
correlation time r c = 1/v, i.e. Z(t) is a symmetric zero- 
mean exponentially correlated noise. The trichotomous 
process is a particular case of the kangaroo process || 
with flatness parameter <p := (Z 4 (t)) / (Z 2 (t)) 2 = l/(2q). 

At great flatness, q — > 0, which is the case addressed in 
the present paper, the transition matrix of our trichoto- 
mous process takes the following form 



0(q). 



(2) 

Within the framework of the three-level noise models 
used in Refs. Hf?j flatness is determined by a parame- 
ter A, which regulates the relative amount of time spent 
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FIG. 2. (a) The surface of current reversals J(D, a, v) — 
for a fixed asymmetry parameter d = 0.0045. (b) The pro- 
jection of the surface onto plane (D,v). The level curves 
correspond to the following values of the noise amplitude: 
(1) a = 27.700000, (2) a = 27.733329, (3) a = 27.753662, 
(4) a = 27.755045, (5) a = 27.765482, (6) a = 27.775920. The 
effect begins at v A = 642.4480 and ends at v c = 654.2591; 
whereas at vb — 646.9504 there occur two two- fold zeros. 

in the state z = as opposed to the states z = a and 
z = —a. If A — ► 0, the flatness ip — > oo and the leading 
order terms in the transitions matrixes of the noise pro- 
cesses of Refs. [^,0 and Eq. (||) become equivalent. Thus, 
at great flatness our trichotomous noise practically coin- 
cides with the noise used by Bier |J and Elston-Doering 

We describe overdamped motion of Brownian particles 
in dimensionless units by the Langevin equation 



dX 
~dt 



h(X)+i(t)+Z(t), h(x) 



dV(x) 
dx 



(3) 



where V(x) = V(x + 1) is a periodic spatial potential 
of period 1. The thermal noise satisfies (£(£)) = and 
(£(£i)£(^2)) = 2D5(ti — ^2)5 where D is the thermal noise 
strength which will below simply be called the temper- 
ature, for the sake of brevity. As said, we take the ran- 
dom force Z(t) to be a zero- mean trichotomous Marko- 
vian stochastic process [PiP Ul^j . To derive an exact for- 
mula for J, the Fokker-Planck master equation corre- 
sponding to Eq. (J3J) is used, supposing that the potential 
V(x) in Eq. (f|) is piecewise linear (sawtooth-like) and 
its asymmetry is determined by a parameter d G (0, 1), 
with d = 1/2 for symmetric V(x). The force caused 
by the potential is h(x) = b := 1/d for x G (0, d) and 
h(x) — —c := —1/(1 — d) for x G (d, 1). Under these 
assumptions, a complex exact formula as a quotient of 
two eleventh-order determinants can be derived for the 
probability current J. To obtain a more manageable for- 
mula, on the assumption that the flatness parameter is 
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FIG. 3. (a) The necessary and (b) the necessary and suffi- 
cient conditions for the four- current-reversal effect. The dot- 
ted region in (a) displays the possible range of the four zeros 
of the current J(D) at different values of d. 



large, cp = l/(2q) ^> 1, we can expand the current in a 
series J — qJ^ -\-q 2 J^ + . . . One of the present authors 
derived for the leading order term q a rather lengthy 
analytic expression, i.e., Eqs. (31)-(33) in Ref. Herein 
we shall apply graphical analysis to study the analytic 
expression of q in the intermediate regimes. 

First we will briefly review the asymptotic limits of the 
current J as the function of v and D, found in Ref. [[)). 
In the case of large flatness at the asymptotic limits of 
both small and large v the function J = J(y) is always 
positive. Hence, in the case under discussion there can 
exist either none or an even number of current reversals 
with v. Two current reversals vs v (as well as D) were 
addressed by Ref. ||. The case of four current reversals 
vs v with the necessary and sufficient conditions for their 
existence was extensively discussed in our earlier paper 

©• 

What concerns current reversals as functions of D, 
then in the asymptotic limit of high temperature, D — > 
oo, we find that 



J: 



q(b — c)a 2 
ISObcD 4 



(4) 



Thus, at high temperatures the behavior of the function 
J = J(D) is uniform: the current is always positive and 
decreases monotonically to zero as D — > oo. At low tem- 
peratures, in the cases a < c and c < a < b the current is 
positive for all values of v and d. For a > b the current 
behaves asymptotically as 



J ~ qv 



D u/c(a-c) 



-v/b(a+b) 



e v/b(a-b) _ e -v/c(a+c) 



(5) 
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FIG. 4. Three current reversals vs temperature D for 
fixed d — 0.005 and a — 255.86. Switching rates 
are (1) z/i = 3.17379365, (2) i/ 2 = 3.17381112, and 
(3) 1/3 = 3.17382860. Curve (2) has three single zeros. Curves 
(1) and (3) have one single and one two-fold zero. 

Notably, J is positive in the case of a < be at any v. If 
a > be, then the current reverses to negative at a critical 
value v = v$. The point of reversal z/q, being a non- 
trivial solution of the transcendental equation J(z/q) = 0, 
can be found by numerical calculation from Eq. (g). In 
the limit of low temperatures, D —> 0, the behavior of 
the function J = J(D) is not uniform. Depending on 
the values of the remaining parameters of the system, at 
D = the function J = J(D) may have a (finite) pos- 
itive, a (finite) negative or a zero value and may start 
to increase, decrease or remain practically constant for a 
while as the temperature D increases [gj. What is im- 
portant about the asymptotic limits in the context of the 
present paper, where we are interested in the behavior 
of the current in the intermediate domains of the sys- 
tem parameters, is the fact that at the limit of large D 
the function J = J(D) is always positive, whereas at the 
limit of small D it can be either negative or positive, and, 
consequently, there can occur any number, odd or even, 
of current reversals with D. 

Before the present paper not more than two current re- 
versals with respect to D had been found for correlation 
ratchets with a simple sawtooth potential. Moreover, in 
a paper |J by one of the authors it was even argued that 
the possible number of current reversals with tempera- 
ture is either zero, one, or two. 

Next we will examine the four-current-reversal effect 
as a function of D (see Fig. 1). Fig. 2(b) exhibits the 
level curves of zero current, J(D,v\d,a) = 0, for fixed 
d = 0.0045 at different fixed values of a = const. The 
level curves may be considered as functions v = v(D) 
(as well as D = D(v)), with d and a being param- 
eters. In Fig. 2(b) the level curves on the left close 
at smaller finite values of v (not shown), whereas both 
branches of the level curves on the right approach zero 
as v grows. Regarding the upper branch on the right, if 
a < d _1 (l — d) _1 , then D becomes zero only at the limit 
v — > 00, whereas if a > d _1 (l — d) _1 , then D is zero at 
a finite v. In view of this, two types of the level curves 
are distinguishable in Fig. 2(b): namely, the connected 
ones (i.e., curves (4)-(6)) and the ones (i.e., curves (1) 
and (2) comprising two components, viz., a closed curve 
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FIG. 5. The necessary conditions for the three-current- 
reversal vs D effect. The curves a = a(d) and v — v(d) 
represent actually narrow regions in the planes (d,a) and 
(d, v). Especially narrow is the bow depicted by a — a(d), 
e.g., at d = 0.001 the noise ranges from a m in = 1280.330 to 
a max 1280.462. 

and a curve with one end open. There is one very special 
level curve (i.e., curve (3)) which intersects itself at the 
saddle point. The four current reversals vs D effect ex- 
ists at a certain fixed value of ^fi X ed = const if and only if 
there exist d and a for which the function v = v(D',d,a) 
has a local maximum ^ max > ^fixed and a local minimum 
^min < ^fixed (see the vertical dashed lines in Fig. 2(b)). 

By gradually varying a and d one can demonstrate 
that the region of existence of the four-current-reversal 
effect as a function of D shrinks to a 4-point C, which 
has the following coordinates: dc ~ 0.0053305, ac ~ 
25.3379050, D c « 0.1244425 and In v c « 6.2977305. The 
four-current-reversal effect is possible if d G (0,dc) 7 ^ £ 
(ac, 00) and v G (yc^oo) (see Fig. 3). The necessary 
conditions for the existence of the four-current-reversal 
effect are shown in Fig. 3(a) by the shaded regions in 
planes (d, a) and (d, z/), and in Fig. 3(b) by the region 
between the ascending outermost curves in plane (u : a), 
while the intensively shadowed narrow wedge-shaped ar- 
eas in Fig. 3(b) fix the values of d, a, and z/, which are nec- 
essary and sufficient for the existence of the four-current- 
reversal effect as a function of D. 

What concerns the three-current-reversal effect as a 
function of temperature D (see Fig. 4), then the effect oc- 
curs only in a very narrow range of the system parameters 
a and v (see Fig. 5). By varying a and d step by step, we 
will see that the region of existence of the three-current- 
reversal effect as a function of D shrinks to a critical 
4-point S in the parameter space, which has the follow- 
ing coordinates: d s ~ 0.03275, D s = 0, a s ~ 38.95953, 
us ~ 3.42833. The three-current-reversal effect is possi- 
ble if d G (0, ds), v G (3.1275,z/5) and a G (05,00) (see 
Fig. 5). 

As demonstrated above, the necessary and sufficient 
conditions for the existence of an odd number of current 
reversals vs D are a > be and v > Vq, where Vq is a 
zero of the function J(u), see Eq. (||). In the opposite 
case, the number of current reversals is even or zero. To 
elucidate the physical meaning of the above conditions, 
let us re-derive Eq. (||) on the grounds of the following 
physical considerations. In the case of great flatness, 
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^ < 1, the noise Z(t) is with overwhelming probabil- 
ity, P s (0) = 1 — 2q « 1, at the state z = and the 
current J can be regarded as the sum J = J + + J - , 
where the positive current J + is caused by the tran- 
sitions z = a <-> z = and the negative current J~ 
by the transitions z = —a <-> z = 0. (Note that the 
transitions z = a <-> z = — a induce current which is 
proportional to q 2 and will be discarded at the present 
approximation.) Under these assumptions the station- 
ary probability distribution at the noise state z = 
with D = consists, evidently, of delta functions at 
x = ±n + d, n = 0, 1, 2, . . . Within the interval (0, 1) the 
center of mass lies at yo = d. Let at the initial moment 
occur the transition z = — » z = — a. The first time when 
the noise turns back to z = is denoted by to- The cen- 
ter of mass stabilizes now at the position y. It is easy to 
find that the center of mass is shifted by Ay (to) = y — yo 
with Ay (t ) = -(n + 1) if (n + l)T~ > t > nT~ + Tf 
and Ay (to) — —n if nT _ + T-f" > to > nT _ . The in- 
terval T~ = Tf + ((1 - d)/(a + c)) is the time that the 
particle for z = —a takes to pass the period L = 1 of 
the potential; Tf = d/(a — b) is the time necessary for 
passing the length d (i.e., from the minimum of the po- 
tential to the maximum in the negative direction). The 
probability W(to) that in a certain time interval (0, to) 
the transitions z = — a — > z = do not occur, is given by 
W(to) = exp (— vto). The probability that such a tran- 
sition will occur within the time interval (to, to + dto) is 
vdto, and consequently, < Ay >= v J °° e~ vt ° Ay (to) dto. 
Considering that the average number of transitions per 
unit time into the z = — ao state is qv, we obtain 
J- = qv < Ay > = -qvW(T^)/[l - W(T~)\. Similarly, 
one can derive the positive component of the current, viz. 
J+= qvW(T+)/[l-W(T+)}, where T+ = (l-d)/(a-c), 
T+ = T x + + d/(a + b). The inequality T+ < T~ being 
equivalent to d < 1/2 , it is evident that the total cur- 
rent J = J + + J~ , whose expression coincides exactly 
with Eq. (^|), can be negative only when Tf < T x + . The 
latter inequality can be written as a > be, which is just 
the necessary condition. 

Following similar trains of thought as above in the case 
of low temperature D <C 1 and replacing the passage 
times T ± and T x by respective mean first passage times 
< T ± > and < >, one can re-derive the necessary 
conditions for the existence of the three-current-reversal 
effect as a function of temperature. Namely, we obtain 
the formulas a = 41/(32d) -0.25 and v = 3.12(1+ 3d) 
which quite well approximate the corresponding curves 
in Fig. 5. 

Finally we will discuss the possible usefulness of the 
phenomena of three and four current reversals. One of 
the first applications of the Brownian ratchet mechanism 
has been to separate particles by size, charge, mass, etc. 
Separation can be achieved even if there are no current 
reversals, as particles of different viscous friction move 
at different speeds. If there is only one current reversal, 
then particles of different friction coefficients, move in 
the opposite directions within the same environment. If 



there occur two current reversals with ^, then particles 
with parameter values within a characteristic interval can 
be separated, as they move in the direction opposite to 
that of the remaining particles. However, it is known 
that the two zeros of the current J(v) generally occur at 
largely displaced values of v. Therefore, in the opposite 
direction move those particles whose friction coefficients 
are within a wide interval and the separation effect is of 
low selectivity. Though this wide interval can be made 
narrower by varying the other system parameters, the 
really narrow intervals occur only in the vicinity of the 
transition regimes (i.e., the transition from two current 
reversals to zero current reversals), where the absolute 
value of the current is small. On the other hand, the 
three-current-reversal with D effect enables us to design 
continuous two-step separation schemes with very high 
selectivity within a narrow v interval. Note that at a 
fixed value of d the range of values of v where the three- 
current-reversal vs D effect exists is very narrow indeed 
(see Fig. 5). 

Within the framework of the calculation scheme of the 
present paper, the absolute value of the net current is 
inversely proportional to the flatness parameter of the 
trichotomous noise, ip = l/(2g), and q as an expan- 
sion parameter is generally considered to be infinitesimal. 
However, we have managed to show by direct numerical 
calculations without using expansion in q that the effects 
are present up to q ~ 0.015. 

Beyond the separation methods, the phenomena of 
three and four current reversals may be of interest in 
biology, e.g., when considering the motion of macro- 
molecules. While it is known that the two-current- 
reversal effect allows one pair of motor-proteins to move 
simultaneously in opposite directions along the micro- 
tubulae inside the eucariotic cells, then existence of three 
and four current reversals will enable such a simultane- 
ous motion of many pairs of motor proteins. Whether the 
intracellular transport makes use of the flashing ratchet 
mechanism as one presently tends to think, or the cor- 
relation ratchets mechanism, or a combination of them 
remains to be seen. 

To summer ize, it is remarkable that the interplay of the 
symmetric three-level and thermal noises in the ratchets 
with a simple sawtooth potential generates such a rich va- 
riety of cooperation effects as up to four current reversals 
with temperature as well as switching rate. The results 
are the more surprising because in analogous model sys- 
tems with the symmetric dichotomous noise the current 
reversals are altogether absent. 
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